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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
TECHNICAL MEMORANDUM 1^ 

EXTREME SPEEDS AND THERMODYNAMIC STATES IN 
SUPERSONIC FLIGHT* 

By Klaus Oswatitsch 



SUMMARY 


The increasing importance of high-speed flow leads to similar prob- 
lems in various fields of research which are summarized in what follows. 
Typical of all cases is the conversion of high kinetic energy into extreme 
thermodynamic states with temperatures of several thousand degrees, fre- 
quently connected with dissociation and ionization of the gas involved. 
There is also a characteristic small sensitivity to the processes dis- 
cussed in the case of gases of low molecular weight (light gases) . 

The penetration of meteors into the atmosphere of the earth at 
astronomical speeds results in temperatures higher than those of the 
surface of the sun. Such temperatures may be produced in shock tubes, 
with light gases used as the driving gas. For supersonic fighters the 
problem of propulsion is less difficult to solve than the problem of 
large heating, on the surface and in the combustion chamber. Finally, 
for the space -travel rocket, astronomical speeds have to be reached which 
require the lightest possible gases as propellants. Here again, dissocia- 
tion. processes in the combustion chamber are of considerable importance. 


1. INTRODUCTION 


If science deals with extreme conditions, it may appear, at first 
glance, as tho ugh the record-crazy present is penetrating into the domain 
of serious research. On the one hand, this is actually the case. Research 
cannot and should not keep aloof from the problems of extreme technical 
questions, and the following expositions are meant to be a critical evalua- 
tion of publicly discussed related pla ns. On the other hand, the study of 

*" Extreme Geschwindigkeiten und thermische Zustande beim 
Uberschallflug." Zeitschrift fur Flugwissenschaften, vol. k, issue 3 
1956, pp. 95-108. Friedr. Vieweg &. Sohn, Braunschweig, Germany, Publisher. 
(Synopsis of the lectures given at the Foreign Institutes of the Technical 
Universities of Vienna and Graz in June 1955 and of the lecture at the 
main meeting of the DVL in Munich on September 30 , 1955*) 
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extreme conditions holds a special lure for research because it leads to 
new and unexpected effects. Regions are reached which, at first, chal- 
lenge engineering and physical intuition horn of experience. Obstacles 
appear in unexpected places while expected ..difficulties often vanish. 

The arrogance which makes us push on into hitherto untouched fields turns 
into humility in the face of the unknown and the unexplored. 

In the following discussions we attempt a maximum of clarity with a 
minimum of mathematical and physical complexity. The results given can 
be found in various reports in various issues of the international 
scientific literature. Only a few papers are quoted explicitly. A 
specialist in this field will gain in what follows - aside from a new 
compilation - at most a few new approximation formulas. However, the 
author hopes to find his most appreciative readers among the experts in 
neighboring fields and among interested novices. Exact derivations have 
frequently been added in small print. “ ' 


2. ASTRONOMICAL SPEEDS AND THERMODYNAMIC STATES 


Table 1 shows, in addition to the velocity of sound c and the 
velocity of propagation of light, approximate values for astronomical 
speeds. In order to have a relationship to the customary velocity scale 
of modern engineering, the velocity of sound is given also in kilometers 
per hour (km/h). 

Compared to technologically realized speeds, we deal in the universe 
actually with "astronomical" numbers. Modern pursuit planes attain sonic 
velocity and will be considerably faster within the next few years. Rifle 
projectiles, antitank and antiaircraft shells as well as naval-gun pro- 
jectiles have initial speeds of about three times the velocity of sound, 
and modern liquid-fuel rockets reach seven times, perhaps even ten times 
sonic velocity. However, when the earth, on its path around the sun, 
encounters a piece of matter, this latter enters the earth's atmosphere 
with a velocity of JO km/sec or about 100 times the speed of sound. A 
shooting star or a meteor restats. The air is not capable of avoiding 
it in time, since a small pressure disturbance is propagated only with 
sonic velocity which is far exceeded by the meteor flying at high super- 
sonic speed. The meteor pushes the air in front of it .in a cushion of 
high density. Ahead of the meteor, density, temperature, and pressure 
Jump, shocklike, in a bow wave (fig. 1) to a multiple of their initial 
values. We speak of a compression shock or, abbreviately, of a shock. 

Here the importance of the ratio of the velocity W and the veloc- 
ity of sound c becomes clear. This ratio is generally designated as 
Mach number M, after the Viennese physicist and philosopher Ernst Mach 
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M = g (1) 

E. Mach was one of the first to recognize the significance of this pro- 
portionality factor for flow problems. 

The pressure increase in the perpendicular part of the bow waVe can 
be easily determined approximately. We visualize ourselves placed on the 
meteor as observers. The pressure increase in the bow wave originates by 
the air, flowing against the meteor with enormous power, being suddenly 
brought to an almost complete stop. The momentum per unit volume is 
pW (p = density) . A momentum stream pW x W flows through a unit area 
per unit time. The rate of loss of this momentum is equal to the static 
pressure p immediately behind the front of the bow wave 

PW 2 = § (2) 

It is true that static pressure exists also ahead of the bow wave, and a 
small momentum stream is present also behind the bow wave since the air 
must flow off around the body. However, both quantities named are much 
smaller than the effects included in (2). 

For ideal gases, there exist the following speed-of -sound formulas ; 
(t = absolute temperature, k = Cp/c v ratio of the specific heats for 
constant pressure Cp and constant volume c v , m = molar weight, 

R = universal gas constant) 

= ^ = . (3) 

Ahead of the shock, the air is certainly to be regarded as an ideal gas. 
Hence there follows from equation (2) with equations (l) and (3) the 
following formula for the pressure rise in the bow wave for high Mach 
number 



(4) 


Formula (^) shows clearly the importance of the Mach number; even though 
it is valid only for higher supersonic speeds - the so-called hypersonic 
speeds - it has the advantage of not containing an assumption on the gas 
state after the compression. This is important, for in the case of the 
meteor flying at Mach number 100 there results (k is the value for air 
ahead of the bow wave: k = 1.4) a pressure increase of 

£ = 1^,000 

p ’ 

thus an extremely high value. 
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An exact calc ulat ion leads to the following formula: Since an equal mass flows in ahdad of the shock 

front as flows off behind the shock front (see, for instance (ref. 1), p. 20 ) , the continuity condition 


pft SS ptf 


(5) 


applies. The exact momentum equation reads 


From equations (5) and ( 6 ) follows 


$ + pft 2 o p + pW 2 


( 6 ) 


<71 

Whereas the last term has do significance whatsoever for high speeds, the parenthesis In equation ( 7 ) In 
the given example leads, according to table 2, to a correction of about 10 percent ccmpared to equation (4). 


In order to form now also an opinion about the temperature increase 
occurring in the how wave, we shall go hack - as in all temperature 
problems - to the energy theorem. Let us consider the unit mass of a 
gas particle which passes through the front of the how wave (fig. 2) . 

The kinetic energy W^/2, and also the internal energy of the unit 
mass e are released; furthermore, the gas flowing out at the pres- 
sure p performs work of the magnitude p/p by displacement of the 
volume l/p per unit mass. On the other side of the bow-wave front 
corresponding energy increases have to he considered; hence there results 
the energy balance 


e 




+ £ + 


W2 

T 


( 8 ) 


This is valid under the assumption made here and always made later, 
that no energy radiates toward the sides; however, at very high tem- 
peratures, this is probably only conditionally true and will probably 
require a correction. In equation (8), the kinetic energy behind the 
bow wave ft 2 /2 is, in practice, quite insignificant. If the energy 
consideration is set up, instead of for the. state behind the bow-wave 
front, for the state at the stagnation point, there even applies exactly 
(with the subscript 0 for the state of rest W = 0; e = eQ, p =■ pq, 

P = P 0 ) 


e o 



e + £ + 
P 


w ® 
2 


( 9 ) 


For an ideal gas there applies, furthermore, the equation of state, 

f'-“ =('T-=v ) 1 * <“» 

and it will he shown that, in what follows, the gases may very well he 
regarded as "ideal.” 

NACA reviewer's note: Although the first equality of equation (10), 

•jj = is applicable in the general case of varying molecular weight, the 

second equality, ~ = (cp - c y )T, is not. 
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Furthermore > if we assume , in order to arrive at numerical values, 
that we are dealing with an ideal gas of constant specific heat, there 
applies for the internal energy, aside from an additive constant, e = c v T. 
Hence it follows, together with equation (10), that 

T j2 

°p T 0 = c p T + \ 

or, introducing the velocity of sound with equation (3) and the Mach num- 
ber with equation (l) 


^ = 1 + -- k M 2 (11) 

T 2 

Hence, for air (k = 1.4) at a meteor Mach number M = 100, it follows 
that the temperature ratio Tq/T = 2,000; for an absolute temperature 
of 225 ° as prevails, for instance, in the stratosphere, this would cor- 
respond to a rest or stagnation-point temperature of Tq = 450,000°. 

At such a temperature, however, the air has certainly been dissociated 
or ionized long before. The assumption of the air as being a gas of 
constant specific heat, which forms the basis of equation (ll), has cer- 
tainly been grossly violated. This is equally true for the state imme- 
diately behind the shock front where the air also has come almost to a 
standstill and where, therefore, almost stagnation conditions exist. 

Whereas the momentum considerations with equation (4) lead, there- 
fore, to a quite serviceable estimate for the pressure increase in the 
bow wave, a more exact knowledge of the gas state is necessary in order 
to arrive at results for the temperature. Figure 3 shows the composi- 
tion of the air as a function of the temperature for normal density 
according to quantum-theory calculations of Burkhardt (ref. 2) .* According 
to these calculations, the composition of the air of the normal state of 
approximately 20 percent O 2 and 80 percent N 2 is maintained up to 
beyond T = 2,000° absolute. Aside from a slight formation of NO, N 2 
and O 2 start dissociating to N and 0 at about 3*000°, and at 10,000° 
the air is a monatomic mixture of nitrogen and oxygen. The mean molar 
weight m has varied inversely as the number of molecules and has at 
10,000°, therefore, only half the magnitude of the normal state. In 
cases of heating beyond 10,000°, the gas is ionized. Due to the forma- 
tion of electrons, the molar weight decreases further. Figure 4 shows 
the internal energy of m = 29g air. This is, therefore, the internal 
energy of a mole at low temperatures and, according to the kinetic gas 
theory for molecules with five degrees of freedom, thus for air of. low 
temperatures , 


me = 5T 

is valid. 

*NACA reviewer’s note: Since the calculations of Burkhardt are now 

outmoded, figures 3 and 4 are only approximately correct. 
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Since the molar weight decreases at higher temperatures ; me is, 
at higher temperatures, no longer the internal energy of a mole. Fig- 
ure 4 shows that the internal energy increases considerably more than 
for the "ideal gas of constant molar weight." This is caused by the 
large energies required to first split the molecules into atoms and then 
to separate electrons from the atoms; such energies are_comparable to 
those of intensive chemical reactions. 

In the entire range considered, the air is to be regarded as an 
"ideal gas" in the sense that the equation of state ( 10^ applies with 
substitution of the particular molar weight m, which is a function of 
the temperature and also of the density. From the standpoint of gas 
kinetics, this means that the volume of the molecules is always small 
compared to the volume of the space in between and that "the forces of 
mutual attraction are insignificant. The dependence of the internal 
energy on the density is only a result of the dependence of the decompo- 
sition of the gas on the density. The potential energy arising from the 
mutual attraction of the molecules, in contrast, does not play any role, 
just as in the case of the ideal gas . “ ‘ T 

For calculating the change of state in the shock wave, we shall U6e" 
a relationship between the thermodynamic state variables which is known 
as "dynamic adiabatic." W and W can be eliminated from the three equa- 
tions (5), (6), and (8) which are valid for all media, and 


§ - e = i(p " § (£ + p ) ( 12 ) 

results. For the extreme states here treated, p » p and, according 
to figure 4, likewise e » e, -whence we obtain the following equation, 
of good approximation: 

£ i = 2S0 . 

p T 

and, with use of the equation of state of ideal gases (10) the approximation 


£ _ 1 _ 2em = 2em fi. _ v 13 \ 

P BT~ if" m 

The difficulty lies in the fact that the state after the compression is 
desired but that there the internal energy e is dependent on T and (5 
in a £omplicated manner. Therefore, we shall simply assume the state pf 
and T and calculate the Mach number of the meteor as well as its flight 
altitude, thus the air density p in which it is flying. 



NACA TM 1434 


7 


A 

The second column of table 2 gives the assumed temperature T 
multiplied by the absolute gas constant. We put there R = 2 cal /degree, 
not 2 cal/g x degree, because the molar weight has been introduced with 
the dimension g (gram) in figure 4. Columns 3 and 4 follow from fig- 
ures 4 and 3, column 5 from formula (13) . It is remarkable here that the 
density in the bow wave increases only by a factor of 6 to 13. This is 
of course not surprising since it is known that the density in the com- 
pression shock of an ideal gas of constant specific heat of k = 1.4, in 
the limiting case of high Mach number, increases only to 6 times the 
initial value. The dissociation and ionization of the gas thu6 causes an 
increased compression and a decreased rise in temperature - as will be 
shown later - and has no important influence on the pressure rise. 

As a result of the slight compression of the gas, an 'ideal gas may 
still be regarded as "ideal" even after the compression shock - however, 
with variable molar weight m, dependent on the state. 

Since j$ was assumed as the normal density, column 5 also determines 
the density ahead of the bow wave and thus the flight altitude of the 
body. It is represented with the assumption of the so-called "standard 
atmosphere" (ref. 3) • In every case, we are dealing with flight altitudes 
in the stratosphere, thus with altitudes which are of interest not only 
for meteors but also for missiles flying at very high speeds. There a 
temperature of approximately T = 223° absolute prevails. Thereby the 
temperature ratio also is fixed. With the equation of state of Ideal 
gases (10) ahead of and behind the front of the bow wave there then fol- 
lows the pressure ratio 


A A . 

p = T p m 

p T p m 


(14) 


and with equation (7) (or with the less accurate equation (4)) the Mach 
number . Because of lack of the m-values for T = 30,000°, the calculation 
in table 2 is limited- to Mach numbers M ^ 44. On the other hand, a 
meteor in the upper stratosphere is strongly decelerated so that table 2 
is of interest for this astronomical phenomenon, too. We shall encounter 
the Mach -n umb er range between 15 and 30 below also in the problem of 
satellites. For this reason, the results are of technical interest as 
well. Of course, they are of importance only for the pressure region 
behind the bow wave. Considerably weaker effects would occur on pro- 
jectiles with conical points. 


As mentioned before, the pressure rise depends in practice only on 
the Mach number. From M = 30 (line before last in table 2) there fol- 
lows with equation (4) p/p = 1,260. In equation (l4) the left side is 
thus fixed with M. The compression Is approximately twice as high as 
for an ideal gas with k = 1.40 and the molar weight only half that in 
the Initial state. Hence the temperature rise in this case is only about 
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l/4 of that calculated with equation (ll); thus, instead of 
'f/T ~ Tq/T = 180, only 45 . For an approach -flow temperature of T = 223°, 
this means, instead of a temperature of T = 40,000°, only 10,000° behind . 
the bow wave - a value which, of course, is still far above the tempera- 
ture of the surface of the sun. The air is practically dissociated com- 
pletely in this case. 

It is true that we make here the assumption that the time intervals during ■which, the gas stays in the 
state considered are sufficient for adjustment of the thermodynamic equilibrium, the dissociation and ioniza- 
tion equilibrium on which figures 5 and k are based. Experiments with the shock tube treated in section k 
have shown that relaxation phenomena do appear. The final state is therefore attained, sometimes, not imme- 
diately behind the shock front but only after passing through a transitional zone. Tt must be taken into con- 
sideration, though, that a delay in the disintegration of the molecules causes a superelevation of the tem- 
perature. The tendency of the gas toward dissociation and ionization is therefore greatly heightened in this 
transitional state which makes the gas there tend toward thermodynamic equilibrium the more rapidly. 


3. DRAG PROBLEM 


The drag of bodies flying at such high speeds is governed exclusively 
by the extraordinarily high pressures in the neighborhood of the stagna- 
tion point. These pressures amount to a multiple of the approach-flow 
pressure. In the face of these facts, even a complete vacuum behind the 
maximum thickness of the body cannot produce a noticeable suction. Fric- 
tion forces can exert a certain amount of influence only in the case of 
very slender bodies which will not be considered here. 

Since the pressure at the stagnation point lies above the pressure £ 
behind the perpendicular part of the how wave (the pressure at the maximum 
body thickness, however, is considerably smaller than £) the drag D will 
he equal to the cross-sectional area F times a mean pressure p = CfJ 
where C is of the order of magnitude of unity. Hence, there follows with 
equation (2) 


D = Fp = CFp = CFpW 2 . (15) 

Thus the drag coefficient cjj referred to the "frontal area" F is 


C-n — 


FpW 2 /2 


- 2C 


(16) 


This result is noteworthy. It indicates that, in the range of very 
high Mach numbers, the drag coefficient becomes independent of M«.. The 
author generally proved this result (ref. 4) for ideal gases of constant- 
specific heat where it was shown that, for blunt bodies, the independence 
of M appears for » 1. The tests (ref. 5) show indeed that the 
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coefficients for a sphere and for a cylinder in an axial approach flow, 
as early as for M = 3, approach an asymptotic final value which lies at 
cp = 0.95 for the sphere and at cp = I .65 for the cylinder (fig. 5 ). 

The corresponding factors C = p/p are C = 0.48 and C = 0.82. Another 
conclusion is surely permissible here. Since the pressure increase in the 
perpendicular part of the shock depends only on the Mach number M and is 
only insignificantly influenced by the dissociation of the gas, the effect 
of dissociation and ionization on the drag coefficient of the gas certainly 
is also slight. This must apply all the more, .the slenderer the body, 
because then the gas also deviates less from the state of constant specific 
heat. The influence of k on cp, too, can only be slight, according to 
equation (l 6 ). 

At this point we want to discuss briefly the problem of drag of a 
meteor where the pure kinetic energy is compared, on the one hand, to the 
work done against drag in penetrating the atmosphere, and on the other 
hand, to the heat of fusion of iron, the chief constituent of most meteors. 
Referred to the unit mass, the kinetic energy is simply W ^/ 2 (table 3). 
With equation (l6), the resistance work in a distance L, at the den- 
sity p, for a sphere of the radius r (a meteor is to be regarded, 
approximately, as such a sphere), is 

DL = Cpr^rtpL ^ = O.^r^pL ^ 

The resistance work for a unit mass then is, with Pg as the density of 
iron 

DL = Q>T1 1 £L W^ 

(4/3)^^ r PE 2 

The height of the "constant-density atmosphere" -(p = 1.293 X 10~3 g/cm.3 ^ s 
8 km. Since the product pL is not dependent on the assumption of p 
(the weight loading per unit area of the ground is the same in any case) we 
may assume the estimate with a mean constant density to he quite usable. 
With Pg = 7*8 g/cm3 there then follows the second value of table 3- 
The radius of the meteor r, expressed in cm, is to be substituted into 
the formula . 

Third in the table , finally , we find the heat of fusion of iron. As 
a quantity having the dimension cal/g, it may just as well be given in the 
form of a kinetic energy per unit mass. Here and later on, heats of trans- 
formation are thus expressed mostly in the dimensions [l/ 2 (km/sec)^| , 
that is, they- correspond to the kinetic energy per unit mass of a body 
when its velocity iB counted in km/sec. Thus 
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jcal/gj = 0.4186 X 10® [(cm/sec)^) 

= 0.4186 x 10“2 |(km/sec) 2 J > 
= 0.09150 2 jJ(km/sec) 2 J 


(IT) 


1 cal/g corresponds, therefore, to the kinetic energy per unit mass of 
a body flying with a velocity of 0.09150 km/sec. This at firBt unusual 
manner of expression will prove advantageous for energy considerations 
in the field of rockets. 


For a meteor of a radius of r = 94 cm or a diameter of 2m, the 
kinetic energy and the resistance work are about the same, that is, the 
boundary where the meteor is brought to a complete stop by the earth's 
atmosphere lies at this diameter. It is true that the pull of gravity 
has not been taken into consideration here, but in the case of a meteor 
it is of little importance. The kinetic energy of a meteor flying in at 
30 km/sec iB approximately 2,000 times (302/0.41 = 2,200) its heat of 
fusion. A spherical body measuring many meters will thus penetrate the 
earth's atmosphere but will not only melt but also evaporate at the 
Impact. A small spherical body generates on its way through the atmos- 
phere - due to its resistance work - so much heat that even a small part 
of this heat is sufficient to melt it. It is therefore not surprising 
that at least the point of maximum temperature rise - the head of the 
meteor - must hielt. 

The considerations on resistance work are associated especially with 
conditions for meteors; thus they do not permit immediate conclusions as 
to the processes in the case of projectiles or rockets. For these latter, 
the drag values per unit mass are sometimes considerably smaller. 


4. SHOCK TUBES 


The considerations following next will be devoted to the problem of 
producing extreme states in the laboratory. For the various customary 
types of supersonic wind tunnels, one encounters many difficulties in the 
attempt to achieve hypersonic velocities. The problem of the energy 
required is not so difficult. Considerable pressure gradients also can 
be achieved; however, the appearance of very high temperatures poses 
technical problems which are very hard to solve. In recent years, a 
test arrangement which, in English, is called "shock tube" has attained 
importance. We deal here with the production of a one-dimensional 
explosion in the laboratory so that the instrument should perhaps best 
be called an "explosion tube" in German. 
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First, let us consider the unsteady propagation of a two-dimensional 
sound wave (fig. 6) . If we deal not with a wave-shaped but with a step- 
shaped disturbance, the sound wave produces a small disturbance of the 
thermodynamic state and of the velocity of sound c and causes a change 
in the flow velocity W. In the reference system of the sound wave, the 
connection between temperature variation and variation of kinetic energy 
was previously indicated in the derivation of equation (ll) . This con- 
nection applies for any stationary disturbances, compressions or expansions 

W dW = Cp dT * 


or, with introduction of the velocity of sound with equation (3) 

W dW § — c dc (18) 

k — 1 

However, we now deal with the disturbance in a sound wave. In order to 
keep the latter stationary in figure 6, the flow velocity must be chosen 

equal to the velocity of sound. Hence, for the change of state in a 

sound wave there applies in the coefficients of equation (18) : W = c. 

The change of state in a sound wave of an ideal gas of constant specific 
heat is given exactly by 

dW § — dc (19) 

k - 1 

Thus, the velocity varies in a sound -wave in the case of air (k = 1.40) 
five times as much as the velocity of sound. For equation (19) it no 
longer matters -whether the sound -wave is considered in a coordinate 
system fixed to the wave step or in an arbitrarily moved coordinate sys- 
tem, whether the wave step travels on or remains stationary; a velocity 
difference dW is independent of the selection of the coordinate system. 
A wave crest, for instance like the one sketched in figure 7> with a 
state of rest (W = 0, c = cq) on the left side, may be regarded as a 
superposition of sound waves, each of which travels in the state created 
by the preceding wave. In figure 7> for instance, the first wave, 
situated farthest to the left, runs with the velocity of sound of the 
state of rest cq. The decrease of the velocity of sound in the first 
"partial wave" - which is only an expression of the drop in pressure and 
density - leads, according to equation ( 19 )> to the gas beginning to 
flow slightly to the right. On this gas flowing to the right, the next 
sound-wave step now moves to the left, etc. Thereby we also obtain 
immediately a picture of the deformation with time of the sound wave. 

The absolute velocity of a sound wave moving to the left is W - c. For 

■* - - 

NACA reviewer's note: This equation should read W dW = -Cp dT. 
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a state of rest, W = 0, it runs with the velocity -cj for W = c, it 
remains stationary. By integration of equation (19) , the flow velocity 
of the gas is readily expressed by c. Taking c = Cq for W = 0 there 
applies 

W = - c) (20) 

Whereas thus the velocity of sound in the wave of finite amplitude 
decreases, the flow velocity increases. With equation (20), there 
applies for the absolute velocity of the sound wave 


W - c = 2 — - (c 0 - c) - c = -c 0 + l±-i(c 0 - c) (21) 

For very strong expansion, that is, for very low temperature and 
a very small value of c, the part of the_.wave situated on the right is 
carried off to the right by the gas flow; that is, the wave flattens 
with time. 

In the consideration of an "explosion" in a tube, this appears self- 
evident. Let us start from an initial state in which on the left side of 
a membrane set into a tube (fig. 8) high pressure exists while the other 
side Is pumped out. Shortly after the membrane is burst, a pressure dis- 
tribution similar to the one sketched and a corresponding sonic -velocity 
distribution exist. The further variation with time is that which was 
discussed with the aid of figure 7» Whereas the pressure on the left 
side is reduced by the progressing expansion wave, the high-pressure gas 
flows, explosively, into the low-pressure side. The attainable maximum 
velocity corresponds to a complete expansion p 0, p 0, T 0, 
c ->0; therefore, according to equation (20) 

W = — - c 
max k - 1 0 



It is true that a pressure gradient sufficient for attaining this 
maximum value is hard to achieve so that we must limit ourselves to con- 
siderably smaller values. The gas expands Isentropically (adiabatically) . 
For this there applies, at k = const 



( 23 ) 


or for the velocity to be attained for a pressure ratio p/pQ with 
equation (20) 
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According to equation (2l) , a high velocity of sound at rest is no 
less important for a high, flow velocity W than a large pressure drop. 
This rest velocity cannot be too greatly increased by heating since a 
doubling of c requires, according to equation (3), a quadrupling of 
the absolute temperature . However, precisely this equation (3) shows 
that gases of smaller molar weight have a considerably higher c. For 

equal temperature, hydrogen (Hg: k = 1.40, m = 2 ) has 1 / 29/2 =3*8 times 

the sonic velocity of air (air: k = 1.40, m = 29 ); this phenomenon is 

explained, from the standpoint of gas kinetics, by the considerably higher 
flight velocity of the molecules of smaller mass. 

If the high-pressure side of the shock tube is filled with Hg and 
the low-pressure Bide with very highly rarefied air, the hydrogen rushes, 
after bursting of the membrane, with great speed into the low-pressure 
side, displacing the air which is compressed in a shock. Due to the 
large pressure differences, the process can be sketched only with a dis- 
torted pressure scale in figure 9* The state of rest of the air is 
designated by the subscript 1 in order to distinguish it from the 
state of rest of the hydrogen. At the surface of contact of the two 
gases, pressure and velocity must be identical. For Hg, the connection 
between W and p is given by equation (2l) . For a moving shock, how- 
ever, it still has to be derived. So as not to lose any clarity, the 
consideration will be limited to the case of strong shocks treated here. 
For the stationary flow there applies equation ( 2 ), with the flow veloc- 
ity behind the shock f? being rather small, compared to W. 

The shock front need only be considered from a coordinate system 
fixed to the approaching gas in order to arrive at the case of a shock 
moving into a gas at rest, as in figure 9» The shock front then has the 
velocity W, and - in the case of a strong shock - the gas behind the 
shock front also has approximately the same velocity. Since the values 
ahead of the shock are written with the subscript 1, and those behind 
the shock without subscript, there applies therefore for the moving 
strong shock 

w 2 _ P_ = 1 K Pi P_ = c 2 _L 2. 

1 p x p x 1 k ± p x 


( 25 ) 
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For Judging the merit of equation (25), we shall briefly derive the exact formula. It la easy to derive, 
from the exact equations (5) and (6), the equation (26) for stationary flow 

- ' <« 

Hence follows 

w - " ■ p - P)(^ - £) _ <*0 

The velocity difference W - W Is the same for all reference systems moved pet^endicularly to the shock - 
front, and is hence equal to the velocity W at the contact surface in figure 8. 771th the designations used 
there, ve have therefore exactly: 

_ (“I 


Equations (24) and (25) may be interpreted as two relationships 
for the unknown state at the contact surface, W and_p/pj_. W can be 

easily eliminated, and we obtain l ~ 

^ (£9) 

which is an equation for the required sonic-velocity ratio of low-pressure 
and high-pressure gas in the initial state when the corresponding pressure 
ratio Pp/Po the pressure ratio in the shock P/Pi are given. The 

two ratios of the specific heats in the initial state, and k, also 

enter into equation (29). It Is assumed that the hlghlpressure gas expands 
like an ideal gas of constant specific heatj the compression shock, how- 
ever, may be connected with dissociation, ionization, or merely with varia- 
tions of the specific heats if it is only sufficiently strong. 


£i = ££L,/H 

c 0 K - 1 Y P 


Since the pressure sharply decreases even in the high-pressure gas, 
a pressure ratio Pq/Pi as large as possible is necessary to guarantee 

sufficient pressure rise in the shock* A pressure of p q = 300 atm on 


the high-pressure side and a vacuum of p^ =* 0*03 atm on the low-pressure 

side are not hard to achieve in the laboratory for small tube dimensions. 
If we provide, furthermore, a pressure ratio in the shock of p/Pj_ = 280, 

which according to table 2 corresponds to a Mach number M = l4*8, and 
assume * k » 1.40, we can calculate from equation (29) 




J 


0.142 = 1 
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If we use H£ as the propellant for the compression of air, the 
required temperature ratio then is, with equation (3) 


Tq fflp Cpg _ 2 7 2 

*1 *1 Ci 2 29- 


5.4 


If, therefore, the air is maintained at room temperature Tj_ = 288°, a 
heating of the hydrogen to . T 0 = 980° (= 600° C) is necessary, to pro- 

duce the conditions of table 2 at M = 14.8. These are not particularly 
large demands on test technique. 

The calculation contains three further inaccuracies which could 
slightly modify the result . A heating of hydrogen leads to a small 
decrease of k, -which could become noticeable in equation ( 29 ) because 
k - 1 enters. Furthermore, according to table 2, the air is compressed 
from 0.03 times to 0.30 times normal density. This leads to a somewhat 
higher dissociation and a somewhat lower temperature but not to any 
essential changes. Finally, the absolute temperature ahead of the shock 
is assumed to be 288°, thus 30 percent higher than in table 2, which 
results in a higher temperature behind the shock and correspondingly 
hlghp-r dissociation. This effect probably overbalances those named 
first . 


Therefore we produce, in the present case, a temperature of more 
t.Vian 5,000° abs in the shock tube, thus approximately the temperature 
of the surface of the sun. In contrast, the temperature in the adja- 
cent H 2 is very low, namely 


J£=i 1/7 

2^ = (iLA * = ( — 280_\ = 0.56 

^cT V p o/ l 10 , 000 / 


or Tgv, = 350° abs 


This reproduces a quite general property of explosions. In chemical 
combustion, as occurs for instance in powder, the temperature does not 
rise beyond 2 , 500 ° to 3,000°, since then, as in the case of hydrogen 
nnfl oxygen, the dissociation becomes so strong that it consumes all 
further available energy. The strong pressure drop in the explosion 
leads to a considerable cooling of the powder gas; the dissociation need 
by no means find the time to recombine. The air, in contrast, shocklike 
compressed by the explosion, is extraordinarily heated up and represents, 
further along, the hot part of the process. 

Hie possibilities for producing high temperatures in the shock tube 
are not exhausted with what we have described above. When the shock, 
in figure 9, finally arrives at the right end of the tube, the thermo- 
dynamic state is further pushed up by reflection. A narrowing of the 
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tube toward the right likewise leads to energy concentration and to 
increases of temperature and pressure. Regarding these facts, or the 
use of other gases,* we refer for instance to the reports of Kantrowitz 
(ref. 6). 

The achievement of such extreme states with relatively small effort 
must of course be paid for with a disadvantage: the transient nature of 

the phenomenon. The phenomenon which is of interest extends only from 
the front of the shock to the surface of contact with the high-pressure 
gas, thus only over a fraction of the entire shock tube, and it races by 
with a speed of (eq. (25) ) 


W = 



3^0 m/sec = 4,800 m/sec 


If we succeed, therefore, in extending the observed state to O.50 m in 
a shock tube of 5 a length, it lasts at one point l/l0,000 sec. This 
is sufficient for photographic recording. Subjective observation is of 
course impossible. - 


5. RAMJET POWER PLANT 


In the face of the unusual conditions encountered at high speeds, 
there comes to mind the question as to what speeds are altogether attain- 
able with airplanes or other flying bodies. We choose as the starting 
point the thrust equation for a free-flying body. As the loss of momen- 
tum leads, on a surface, to a pressure rise, to a force. in the flow 
direction and, in the case of a body, in the last analysis to a drag 
(shown for instance in eq. ( 15 ) ) , so the production of momentum leads to 
a thrust S. In what follows, we shall consider the simplest and most 
important case in practice: that the production of momentum occurs in 

a region of constant pressure (fig. 10) . It is unimportant whether pres- 
sure differences appear on the thrust body itself. Let the pressure be 
constant only in a certain region characterized by the rectangular bound- 
ary of the body. Let us consider there the entering jet with the cross 
section f ro > the velocity W^, and the quantity G,* flowing through per 
unit time, and the exhaust jet with the corresponding quantities without 
subscript. 

The loss of momentum in a certain direction, either by stopping or 
by rectangular deflection, produces a force K = pWf X W = GW; the pro- 
duction of momentum generates the corresponding thrust. ' Thus there 
results for the total thrust 


S = GW - G ro W ro 


(30) 
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In many cases a definite inflowing mass is used for increasing the veloc- 
ity; this applies exactly for all propellers and with very good approxima- 
tion for all Jet power plants. Then G is equal to the quantity G m 
coming from the infinite approach-flow region and 

S = Goo(W - Woo) (31) 

is valid. 

In the case of the rocket, however, which will he treated in the 
next section, the momentum, is imparted to a mass which is carried along, . 
G^ = 0, and 

S = GW (32) 


is valid. 

The ramjet power plant operates as follows: A part G of the 

approaching air is given an increased pressure hy being stopped. In 
this state of rest the air is heated whereby its density greatly decreases. 
The heated air finally is expanded again to the initial pressure. The 
velocity W attained by the air in the pressure gradient is larger than 
the velocity W ro lost by the air during the pressure rise of the stop- 
ping because the air density during the expansion is smaller. For the 
lighter the air, the greater is its acceleration in a given pressure 
gradient . 

This type of propulsion was discovered by the Frenchman Lorin at a 
time when flying was still done at low subsonic speeds. For supersonic 
velocities, this type of propulsion was rediscovered by Trommsdorff in 
the early years of the last war. As will be shown directly, this type 
of propulsion is of eminent importance precisely for supersonic speeds. 

In order to understand this, one has only to start from the principle 
that the efficiency tj c of an ideal thermodynamic engine is given, 

according to Carnot, by the ratio of temperature rise between heat input 
and heat exhaust to the temperature of the heat input. The heat input 
takes place at the stagnation temperature Tq. Here as well as for the 
gasoline engine, the cooling of the working gas is replaced by an exchange 
of the hot air blown off in favor of newly supplied cold fresh air. With 
the equation (ll), applied to the approach-flow state T^, M^,, there 

follows thus for Carnot's efficiency 
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The approximation of an ideal gas of constant specific heat, assumed 
for equation (11), is sufficiently satisfied for the following considera- 
tions. Table 4 gives a few results. 

Whereas, therefore, the efficiency of the ideal ramjet power plant 
is still low for flights at the speed of sound, it equals, for M^ = 3, 
a diesel engine working ideally with a pressure of 37 atm! With increasing 
flight Mach number Moo "the efficiency improves still further. Thereby 
this type of propulsion proves to be the ideal engine for steady super- 
sonic flight. Of course, same losses are to be expected for this type of 
propulsion, too; among them, the kinetic energy being lost with the thrust 
jet is by far the most important. Thus only part of the performance given 
by Carnot's efficiency is used "profitably." To an observer standing on 
the ground and considering the body flying at Woo , the work performed per 
unit time against the drag forces or against gravity is 


WJ3 = GWoo(W - W„) (34) 

The exhaust jet, however, discharges an amount of kinetic energy per unit 
time given by 


„ (w - w „) 2 

° 2 

The discharged heat energy need no longer be taken into consideration. 
This energy has already been represented by Carnot T s efficiency. The 
ratio of the desired mechanical output and the sum of the two actually 
achieved mechanical outputs yields the jet efficiency 



^st 


2W 00 (W - Woo) = 2 Wqq 

2W 00 (W - Woo) + (W - Wo,) 2 W + W ‘ 


(36) 


The product of the two efficiencies tj s ^. and tj c then gives rather 

accurately the efficiency which has to be taken into account in practice. 
An example will clarify the possibilities. 


With assumption of isentropic (adiabatic) expansion it is easy to 
derive from equation ( 11) the generalized Bernoulli equation 


W 2 = 


r t i 

/ 

K-l 
•n \ K 

1 - - 2c p t o 


SL) 

?o) 

L OJ 




(37) 


which is similar to equation (24), yet essentially different from it. 

At a certain pressure ratio p/pQ, W is proportional to which is 
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a .square-root relationship between W, c, and 1 appearing in different 
forms. If in the case M^, = 3 the stagnation temperature is quadrupled 
by heating, which leads to a temperature of 2,536° abs, the exhaust veloc- 
ity is doubled: W = 2^. We then arrive at efficiencies of 


^st = 2 / 5 11 = Vst = 0,43 

thus at very promising values. Of course, the jet efficiency could be 
improved by less heating. However, the thermodynamic efficiency alone 
is not the governing factor for the designer. Above all, a certain 
thrust S is necessary. 

If the problem consists, for instance, in flying at constant speed, 
in steady flight, for a long distance, say across the Atlantic Ocean, 
the thrust S is equal to the drag D or with equations (31) and (l6) 


GjW . Woo) = fcoPcoWjW - W„) = CjFpJI 2 J2 


Thereby the ratio of the cross-sectional area of the oncoming jet f ro 
and the frontal area of the body F is 

^*00 __ C D Woo 
"F “ 2~ W - Woo 

The drag coefficient cp of a flying vehicle lies far below that 
of blunt bodies (fig. 5)- If if is assumed to be cp> = 0.30 (referred 

to the frontal area, not, as in the case of airplanes, to the wing area 
or the like), 15 percent of the frontal area is required for the oncoming 
jet at W = 2^00$ this appears reasonable. 

It Is true that a reduction of the difference W - W OT would lead 
to an increase of but at the same time it would increase the jet 

area f^ required according to equation ( 38 ); which is undesirable. 

The designer therefore demands a certain "thrust concentration" and fore- 
goes an excessive increase of the jet efficiency. 

It Is one of the email Inaccuracies of the present representation that we always speak of heating but 
that, in practice, this heating takes place by combustion of a hydrocarbon. An input of mass occurs which is, 
however, rather insignificant. For instance, the combustion of g a soline and benzol in air under constant 
pressure in the mass ratio 1:100 yields a temperature rise of approximately ^00°. As for the rest, a pert of 
the smellier neglected effects acts favorably. 
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With the quadrupling of the stagnation temperature at M w = 3 the 
boundary is reached at which the dissociation of the air sets in. There 
would not be any objection against the storing of energy in the form of 
dissociation if this energy would become free again during the cooling 
off of the jet in the expansion. This would even have the advantage that 
the combustion-chamber temperature would not rise above a certain limit. 
The control of high combustion-chamber temperatures represents, anyway, 
one of the serious technical problems. However, the time available 
during the expansion is too short for recombination of the dissociation 
and adjustment to the cooler thermodynamic equilibrium state since a 
stronger cooling occurs only at very high speeds. Thereby a considerable 
part of the energy going into dissociation is lost. 

The temperature up to which heating may be Continued in the combus- 
tion chamber is thus limited at least by the onset of dissociation - 
unless the strength of the combustion-chamber walls (essentially a cooling 
problem) imposes lower limits. This means a barrier for the heat supply 
which restricts application of the ramjet power plant by a limiting Mach 
number. Let us assume the temperature of 2,536° abs - an enormous heat - 
to be the upper temperature limit in the combustion chamber. Then there 

results for = 5 an exhaust velocity of W = 72,536/1,338 W^ = 
compared to W = 2W m for = 3. This means a considerable increase of 
the area ratio f^/F , even though a small reduction of cj) with may 

be assumed. 

However, temperature problems do not occur merely in the combustion 
chamber but on the entire body of the flying vehicle . It is shown that 
the stagnation temperatures given by equation (ll) and calculated in 
table 4 are attained not only by ordinary storing but also - almost - by 
the stopping of the air due to friction at the surface of the flying body. 
Even though these temperatures lie far below the heating temperatures in 
the combustion chamber, the related problems are not any smaller since 
enormous difficulties are encountered in the cooling of larger parts of 
a flying body. If the vehicle is manned, a very radical cooling of the 
cockpit is needed already for M ro = 3* The electronic components of 
remotely controlled bodies also are very sensitive to temperature; they 
must be cooled or strongly heat -insulated. The simplest- solution depends 
essentially on the length of the flying time. In the case of wing con- 
structions, the strength properties of the metals must be exploited to 
the utmost. The metals greatly vary in strength with the temperature. 
Al umi num melts at 660° C. Due to its higher melting point of 1,550° C, 
soft iron is harder than steel at high temperatures . In finer construc- 
tions, however, not only the general strength must be considered. The 
heat transfer to the wing and the heat conduction in the wing may lead 
to a nonuniform temperature distribution and to a warping of the wing. 
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Still another difficulty must "be pointed out which lies in the 
appearance of high centrifugal acceleration. If the latter is designated 
by g and the radius of curvature of the path by R, the following rela- 
tionship, which is known from, mechanics, applies: 


gR = W2 (39) 

For a speed of W = 1,000 m/sec, thus approximately M = 3> and a path 

radius of R = 100 km = 105 m, one arrives at a centrifugal acceleration 
of 10 m/sec^, thus approximately at the acceleration of gravity. If 4g 
is regarded as the maximum admissible acceleration for the human body, 
only path radii larger than 25 km are admissible for M = 3« The pilot 
must turn in a circle of this radius if he misses his target. For the 
wings and control components of a pursuit rocket, the intended path radii 
thus are sometimes much more important than the production of lift for 
overcoming gravity. 

All this shows that the main difficulty in supersonic flight does 
by no means lie in the production of the required forward thrust, but 
in the control of the temperatures and centrifugal forces. 


6 . ROCKET 


In a well-known and very good book for young readers , written in 
1931 * on various philosophical , mathematical, and technical problems 
there Is the sentence: 11 It can easily be calculated that the upper 
limit to which present fuels can take us, is about 400 km above the sur- 
face of the earth • . The corresponding calculation is not given; 
but, doubtlessly, itjLs based on the relationship that the energy required 
for lifting a mass M (M for the mass in contrast to M for the Mach 
number) to the height H can at most be equal to the reaction energy or 
combustion energy which this mass may contain. At the anticipated height 
of 400 km - compared to the size of the earth* s radius of 6,380 km - the 
gravitational acceleration g may be regarded as constant. With q as 
the reaction energy per unit mass there results, in this manner, the 
energy equation: force Mg times distance H equals mass times reac- 

tion heat q. The mass may be cancelled from this equation with the 
result 


' gH = q. (40) 

A good powder has about q_ = 1,000 cal/g or, according to equation ( 17 ) , 
q. = 186 (km/sec)2. With g = 10 m/sec^ = 10**2 km/sec2, equation (40) 

leads to the height of rise H = 4l8.6 km. 



22 


NACA TM lb$k 


As to the reaction energy, one could of course raise the objection 
that much greater heights could be attained with liquid'- fuel rockets (see 
for instance the oxygen-methane reaction in table 6) and that, more 
recently, atomic energies open up completely different possibilities. 
However, we shall not enter into details about this since the expression 
(40) itself is based on a wrong assumption: that the reaction energy of 

a mass is consumed in lifting or else for accelerating one and the same 
mass. The processes have been sufficiently discussed in the preceding 
section to understand them correctly in what follows. 

The thrust of a rocket is given by equation (32) ; W is the "exhaust 
velocity," the velocity in the rocket jet relative to the rocket body. 
After all the propellant has been ejected, the energy balance is such 
that the thermal energy of the propellant has been used to inpart to the 
rocket body a certain kinetic or potential energy. Undesirably, mechani- 
cal energy (kinetic or potential) has, in addition, been given to the pro- 
pellant and possibly to masses released or left behind in flight. The 
distribution of the mechanical energy between the masses involved is 
rather nonuniform and may, theoretically, be completely in favor of the 
rocket . 

Let us assume, for example, that the acceleration of the rocket 
occurs horizontally - or also vertically but at sufficiently slight dif- 
ferences in height - so that the variations of the potential energy of 
the masses under consideration do not play any part. Lit us assume, 
furthermore, that the exhaust velocity W of the propellant jet can be 
regulated and is, in each case, adjusted so that it is equal to the flight 
velocity of the rocket; then, the ejected propellant mass is at rest, 
relative to the observer. In this ideal case the entire thermal energy 
of the propellant mass flj is used for imparting a kinetic energy to the 

rocket end mass Mr. 

It is true that the "thermal energy" there is not identical with the 
"reaction energy" at a particular temperature, as is usually indicated 
in tables (cf. table 6). The variations in the internal energy of the 
propellant also play a certain role. For the considerations of this sec- 
tion, however, which have only the purpose of orientation, these varia- 
tions are unimportant . But if compressed air is used as the "rocket 
propellant," the "thermal energy" consists exclusively in variations of 
"internal energy." If thus the entire theimal energy per unit mass of 
the propellant which has been released is denoted by q and V is the 
flight velocity of the rocket, there applies in the ideal case 

( 41 ) 


y2 _ 

— Mr = q% 
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By selecting a sufficiently small mass ratio we can there- 

fore impart to the rocket an energy far exceeding equation (40). Evidently 
equation (4l) defines an unattainable ideal case whose value consists, 
above all, in establishing a limit for the most favorable conditions. 

Strictly speaking, the exhaust velocity W must not be set equal to the flight velocity, at least at 
the beginning. At the start a W, even though a minimum W, must exist so that a starting thrust may he 
present* Thus we deal in this respect, as well, with a theoretically constructed limiting case. 

In spite of all the defects inherent in the ideal case investigated, 
it teaches clearly a few fundamental principles. Extreme end velocities 
require large masses of propellant. Bnpty fuel containers will have to 
be dropped, wherever possible, in order not to accelerate unnecessary 
dead weight. The most favorable exhaust velocity W lies below the end 
velocity of the rocket. Thus the rocket jet has, at the start, absolute 
velocities which are directed against the flight direction. During the 
flight, the absolute velocities of the jet decrease more and more, van- 
ish, and are directed in the direction of flight at the end of the com- 
bustion time. In this manner a minimum of mechanical energy may be lost 
to the propellant. 

The rocket combustion chamber of constant exhaust velocity W repre- 
sents the case which is by far the most important in practice. We shall 
compare this frequently treated case (ref. 7) with the ideal case. Let 
U and $ be the velocity and mass of the rocket during the combustion. 

The thrust then imparts to the mass an acceleration dU/dt and there 
applies according to Newton 


M = GW (42) 

dt 

an equation which is in terms of velocities. The propellant mass ejected 
per unit time equals the decrease of the rocket mass 

Eliminating dt from equations (42) and (43), we can relate the 
variation . dM directly to the variation dU 

dXJ = dfl 

W " a 

For W = Constant, this can be easily integrated 

- 2 = in m + Constant 



24 


NACA ®4 14^4 


and since for the beginning and the end of the acceleration 

U = 0 l A = + M/p , 

U = V: A = % 

is valid, there follows finally 

V _ in + for or % + % c e v/w (44) 

w % % 

The fraction (Mr + that is, the starting mass of the rocket 

divided by its final mass, is denoted as mass ratio. According to equa- 
tion (44), this ratio must assume the value e = 2.718 (in e =* 1) if 
the rocket is to attain the exhaust velocity W of the propellant, 
unimpeded by air drag or gravity. For determination of the velocity with 

the smallest jet losses, the kinetic energy of the rocket would 

have to be related to the kinetic energy of rocket and propellant in the 
final stage. Thus we must sum over the kinetic energy of the various 
parts of the jet. Without demonstrating here this comparatively elemen- 
tary calculation, we shall only report that the highest value of the thus 
defined jet efficiency results for a mass ratio of 5- We obtain 

approximately 


( 45 ) 


In the second line, the values for V = W are given. Thus the 
ideal case, given by equation (4l), is not so very extraordinary, after 
all. With the usual rocket of constant exhaust velocity W and ejec- 
tion mass G, more than 50 percent of the ideal case can be attained. 

The most serious difficulty lies rather in the achievement of the required 
mass ratio. The construction of a rocket, the propellant mass of which 
amounts to 63 percent (v/w = 1.0) or even 80 percent (v/w = 1.6), repre- 
sents an enormous design problem. For the V2, the approximate value 
(% + %)/ fl R = 3*5 is given in reference 8. Only the smallest part of 

the end mass of the rocket, however, consists of the useful load 

to be transported; the largest part is occupied by the mass of the fuel 
tanks, of the rocket-combustion chamber, and of the various necessary 
accessories. 


(Mr + %)/% - 5 

v/w = 1.6 

%t = °* 6 5 

(% + %)/m t = 2.72 

v/w = 1.0 

list “ °*58 
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For judging the difficulties of building interplanetary rockets, 
we shall first calculate the necessary starting velocities. We assume 
that the body is directly Bhot away from the earth’s surface into a 
vacuum. Actually, the rocket is driven comparatively slowly through 
the lower atmosphere, the troposphere, in order to keep the losses by 
air drag low. Outside the troposphere, the rocket is then rapidly accel- 
erated to the final velocity V. However, since not only the troposphere, 
with approximately 10 km height as a region, of possible friction losses, 
but also the height required to attain V is small compared to the dis- 
tances considered, we may calculate as if the body were entering the 
vacuum with V directly from the earth’s surface . 

The gravitational acceleration g decreases quadratically with the 
distance R from the center of the earth. On the surface of the earth, 
it is R = Rq = 6,580 km, g q = 9*81 m/sec2. Thus 

g = So@) W 

is valid. If we now assume that a body rotates as- a satellite, at the 
distance R from the center of the earth, with the velocity V s , the 
relationship (59) between path radius, centrifugal acceleration, and 
velocity is valid, so that 


V 


s 


2 



(Vf) 


Table 5 shows numerical values and the Mach numbers referred to a 
sonic velocity of c = 330 m/sec. At such speeds a body does, there- 
fore, not need wings. The path curvature produces the required "lift." 

If a body is to fly at the distance R with V s , its starting 
velocity V must be so large that the decrease of the kinetic energy is 
equal to the work performed against the gravitational pull. All these 
quantities are proportional to M; for this reason, the mass is cancelled 
from the calculation: 



T" “ s o R o ( 1 " 2 Ip) (48) 


or with equation (Vf) 
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The kinetic energy required for escaping the fie.ld of gravity of 
the earth (R -» ») is therefore twice the kinetic energy of a satellite 
of equal weight at a small distance from the surface of the earth. For 
the former, we obtain with equation (17) 

= = 6 3(Wsec) 2 = 15,000 cal/g 

This quantity of energy is not unattainable. The oxygen-methane reaction, 
for instance, yields approximately l/6 of it (table 6) . Thus energies 
are concerned which are perfectly producible. However, the velocities 
to be reached are of astronomical dimensions as is shown by a comparison 
of tables 1 and 5 - a hint that the movements of the stars must be traced 
back to thermal energy sources in atomic reactions, not to chemical 
reactions . 

Table 6 shows several reaction heats and the mass ratio for the ideal 
case, equation (4-1). The melting of a body does not change anything in 
its molecular composition and relatively little in its structure; not even 
the density changes essentially. It is therefore understandable that the 
"reaction heat" connected with melting is generally much smaller than that 
connected with a chemical reaction. The oxygen-hydrocarbon reaction 
02-methane, the detonating-gas reaction, and the combination of atomic to 
molecular hydrogen are given as a representative for the last-named reac- 
tion. It corresponds, of course, only to a thermo-chemical equation, not 
to a practical process, since the atonic hydrogen is not stable under 
"normal" thermal conditions. Nevertheless, this hypothetical process is 
very useful for the following consideration. 

In regard to the structure of matter, chemical reactions are changes 
in the shell of the electrons. In the case of light elements and light 
molecules where the electron shell surrounds only few protons and neutrons, 
such changes may require the same energy as in the case of heavy elements. 
Thus, very much higher reaction heats are possible when substances of low 
molar weight participate than in the case of substances of high molar 
weight m. 

Finally we have the heat quantity of an atomic reaction of 6 = 0.005 
mass defect which, of course, amounts to an extraordinarily strong atomic 
reaction. The heat set free is equal. to the mass defect times the square 
of the speed of light. The velocity corresponding to the equivalent 
kinetic energy is thus equal to the speed of light multiplied by \^5. 
Changing the structure of the atoms leads to energy transformations which 
are incomparably larger than those of chemical processes. Of course, for 
atomic reactions much more than for chemical reactions, it is always true 
that only part of the propellant represents really active reacting sub- 
stance while another part is carried along as "dead substance." Table 6 
is not supposed to reflect more than a picture of the possibilities. 
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It becomes clear that even the mass ratios of 6.9 and 5*7 calcu- 
lated for the ideal case (eq. (h-l) ) , required for surmounting the gravity 
pull, are extraordinarily high for a hydrocarbon and a detonating-gas 
reaction and pose probably unsolvable problems to the designer. (As was 
emphasized before, the reaction heat in table 6 and q in equation (h-l) 
are only approximately equal.) The pertaining energies expressed in 
velocity squares show that these values also lie throughout considerably 
below the required values of V given in table 5 • 

It is instructive to look somewhat more closely at the conditions, 
for instance, in the case of a hydrocarbon reaction. From the known 
relationships between the specific heats 


Cp “ Cy = R/m 


“P/Cv = 


= K 


if follows readily that 


r _ k R 
C P “ k - 1 m 


(50) 


Since in the 0 2 , CH^ reaction the number of molecules is maintained, 

the molar weight is not changed, either, by the process and lies, with 
m = 27, close to the value of air. With k = 1.30, equation ( 50 ) leads 
with R = 2 cal/(g degrees) to a value of approximately 
Cp = O .32 cal/(g degrees). Thus the reaction would result in a tem- 
perature increase of 

AT 1 — ft — 2,3^0 _ Q rj(X)0 

^ " 0.32 ' ’ 

On the one hand, the temperature would be completely intolerable; on the 
other, it is never reached, due to dissociation setting in, as shown by 
the example for air calculated at the beginning in connection with the 
flight of a meteor. These enormous energy transformations in the rocket 
combustion chamber are therefore not even desirable and the carrying along 
of "dead substance" is quite welcome. The maximum velocity attainable in 
the jet can be easily given from equation ( 37 ) under the assumption of 
pressure drop to vacuum p -» 0. With equation ( 3 ) , W^y is 

w max = ^ 2 CpT Q = ^ c 0 


(51) 
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The maximum velocity attainable in the steady pressure gradient, e^resaed by equation (51) , la thus con- 
siderably manwr than that of an unsteady two-dimensional explosion, equation (22). There exists, also, an 
essential difference between the two processes. Whereas in the case of the steady phenomenon every particle 
is equal and gradually assumes in the pressure gradient p ->0, the velocity of equation ( 51 )# wily a sma J l 
part of the matter involved is distinguished by attaining the velocity of equation (22) for p -» 0 in the case 
of the two-dimensional explosion. Furthermore, it riruat be noted that the flow velocity in the unsteady two- 
dimensional explosion exceeds the flow velocity of the steady flow only in the case of a high pressure gradient. 
The pressure gradient for which both velocities are equal can be easily determined by equating V/cq according 
to equation (24) and according to equation (57)* We obtain ... 


or for k « 1.40 




Pg 

P 


17.1 


The speed of sound depends essentially on the molar weight and thus, 
for the mixture investigated, does not differ very much from the value 
for air at the same temperature. If we assume, in consideration of mate- 
rial and dissociation, an absolute temperature of 3*000° abs, we have 

approximately c q = l/lO x 350 m/sec = 1.1 km/sec and according to equa- 
tion ( 51 )* with k = 1.3* the exhaust velocity W = » 2.8 km/sec. 

This value, however, lies far below the ideal Value. For leaving the 
region of the pull of gravity, according to table 5 and equation (45)* 

W = 0.62 x 11.2 km/sec = 7 km/sec would be desirable. Somewhat more 
than W = 2 km/sec was attained with the V2 during the first years. 


Hence, it follows that substances of low molar weight are highly 
desirable not only as reacting masses but also as dead rocket-jet masses 
for extreme speeds of travel. Due to their high sonic velocity, they 
also make high exhaust velocities possible, according to equation ( 51 ) • 

The transformation of higher energy quantities in the combustion chamber* 
in contrast, does by no means lead to higher temperatures, compared to 
heavier substances, because, according to equation ( 50 ), the specific 
heat is much higher for smaller m. Neither is the tendency toward dis- 
sociation at all stronger for gases of small molar weight (Cf .* for 
instance* H 2 according to reference 9*) 

The application of an atomic jet, in spite of the enormous energy 
transformations connected with it, does not immediately lead to a solu- 
tion of the difficulties: Exhaust velocities far exceeding the flight 

velocity mean an enormous waste of energy. However, use of atomic energy 
as an energy source for the heating of a dead substance. of low molar weight 
is to be regarded as a solution promising progress. This dead substance 
should, of course* be carried along as a fluid or solid matter in order 
to keep the container volumes tolerable. For this reason it is very 
regrettable that H 2 , as the lightest substance by far, has a boiling 

point of - 253 ° C and* as a fluid* in addition, has very low density. 
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A possibility of achieving the necessary mass ratios for a flight 
into interstellar space or for sending out an artificial earth satellite 
consists in shooting a "daughter rocket" from a "mother rocket" and pos- 
sibly repeating this procedure in several stages. It is evident that 
the greatest part of the propellant energy utilized must be used to 
accelerate the entire vehicle and only a small part can be used to the 
advantage of the basic purpose. Thus, it will have to be regarded as 
an extraordinary technical performance if we succeed, within the next 
few years, in shooting a satellite of 50 kg weight into an altitude of 
several hundred km. Within a foreseeable time, the shooting of space 
rockets will surely not become such a familiar phenomenon as the pene- 
trating of meteors into the earth’s atmosphere treated at the beginning. 
Certainly it will always be a very costly procedure; it is surely extremely 
optimistic to insert for the actual pay load an expenditure of 10 times the 
amount calculated In connection with equation (48), thus 150,000 cal/'g. 
With this digression to the questions of interplanetary aviation we shall 
close the section on rockets. 


7. RESUME 


All problems treated have in common the high-temperature differences 
connected with high supersonic speeds - a consequence of the large kinetic 
energies inherent to such flows. High temperatures appear, therefore, 
when the air in front of meteors and on fast-flying bodies is stopped. 

High temperatures are necessary when high speeds are to be produced in 
propulsive jets of flying bodies. In increasing the temperature beyond 
several thousand degrees, the heating is greatly reduced by dissociation 
and ionization phenomena since a considerable part of the energy is used 
for splitting up the molecules Instead of for increasing the temperature. 
The influence of the molar weight is very great . The light gases, like 
hydrogen, with the considerably higher flight velocities of the molecules , 
permit the achievement of high flow velocities more readily than gases of 
high molar weight. With hydrogen as the propellant it is therefore possi- 
ble to produce, in the laboratory, transient high temperatures in air, and 
gases of low molar weight appear particularly suitable as the propellant- 
jet substance for rockets. With them, the astronomical velocities which 
are necessary for space rockets are most nearly attainable. 

The high centrifugal forces, which balance the gravity pull at the 
high velocities in space, make themselves clearly felt even at supersonic 
speeds. They impose a limitation on the path curvatures of flying bodies. 


Translated by Mary L. Mahler 
National Advisory Committee 
for Aeronautics 
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TABLE 1.- ASTRONOMICAL SPEEDS 



km/h 

km/sec 

Velocity of sound c in the stratosphere 
(temperature -50° C) 

1,080 

0.30 

Circumferential velocity of the earth at 
approximately 50° latitude 

1,080 

0.30 

Velocity of revolution of the moon about 
the earth 


1 

Velocity of revolution of the earth about 
the sun 


30 

Velocity of the sun with respect to the 
next constellation 


20 

Velocity at the edge of the universe 


50,000 

Velocity of light 


300,000 


TABLE 2.- COMPRESSION-SHOCK WAVE FOR 0 = NORMAL DENSITY 


A 

T, 

abs 

RT, 

&al] 

A 

em, 

[cal] 

m/m 

P/P 

Flight 

altitude, 

[km] 

T/T 

p/p 

M 

(for 

T = 

223°) 

1 x 103 

2 x 103 

5 x 103 

1 

6.0 

14.3 

4.5 

27 

4.8 

5 X 103 

10 X 103 

56 x io3 

1.20 

10.4 

17.8 

22.4 

280 

14.8 

10 x 103 

20 X 103 

235 x lo3 

1.98 

12.9 

19.2 

44.8 

1140 

29.8 

20 x io3 

40 x 103 

500 x lo3 

2.48 

11.1 

18.2 

89.6 

2460 

44.1 
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TABLE 3.- COMPARISON OF THE ENERGY FOR VERTICAL 
PENETRATION OF THE ATMOSPHERE 


Kinetic energy per unit mass: 


Resistance work per unit mass: 


. V/2 


9b W 2 


r [cm] 2 

Heat of fusion of Iron: 49 jcal/gj = 0.4lj^i(km/sec)^j 


TABLE 4.- CARNOT EFFICIENCY t^, CORRESPONDING PRESSURE RATIO OF 
THE PISTON ENGINE, AND STAGNATION TEMPERATURE 
OF THE RAMJET POWER PLANT 


Moo 

n c 

p'o/p 

ToAoo 

T«, = 22 3 

0 abs 

T 0 

0 * 

1 

ro 

a 

1 ' 

1/6 

1.9 

1.20 

268 ° abs 

-5° c 

3 

9/14 

37 

2.80 

634° abs 

361 ° C 

5 

5/6 

530 

6.00 

1,338° abs 

1,065° c 
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TABLE 5.- STARTING AND SATELLITE VELOCITIES FOR VARIOUS 
DISTANCES FROM THE EARTH f S SURFACE 


R - Rq, 
[km] 

[km/sec] 

v, 

{km/sec] 

M 

for 

c = 330 m/sec 

0 

7.9 

7.9 

24 

500 

7.3 

8.2 

25 

1,000 

6.8 

8.5 

25.5 

2,000 

6.0 

8.9 

27 

CO 

0 

11.2 

34 


TABLE 6 .- REACTION HEATS 



[cal/gj 

|^(km/sec) 2 J 

Mr + Mrp 
®R 

Heat of fusion of iron 

b9 

0.64 2 


Black powder 

630 

2.3 2 

25 

O 2 and CHij. 

2,560 

4.6 2 

6.9 

Og and H 2 -fr^O 

3,220 

5-2 2 

5.7 

2E -»H 2 

l /2 percent mass defect 

49,200 

20 . 2 2 
30 , 000 2 

1.3 











Bow wave 


Figure 1. - 


Figure 2. - 



Bow wave for -high. Mach number. 



3ow wave in the stagnation region. 
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Figure 3. - Composition of the air for normal density, as a function 
of the temperature (according to G. Burkhardt (ref. 2)). 
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Figure 4. - Internal energy of m = 29g air for normal density, as a 
function of the temperature (according to G. Burkhardt (ref. 2)). 
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Figure 5. - Experimental drag coefficients for a sphere - and a cylinder 
in axial flow (according to ref. 5) . 

Peenemunde measurements 

Gdttingen measurements 


• • 
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Figure 8. - Initial state in the shock tube. 1 



Figure 9.- Explosion process in the tube. 
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figure 10. . m, 

thrust body. 


Mc< ■ Un «'*y Fi^d, v,. 


